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I obtain the inverse of the correlation length exponent at 
the superfluid-Bose glass quantum critical point as a series in 
small parameter yd—l, with d being the dimensionality of 
the system, and compute the first two terms. For d = 2 I find 
u s = 0.81 and v c = 1.03, for short-range and Coulomb inter- 
actions between bosons, respectively. When combined with 
the exact values of the dynamical critical exponents, these 
results are in quantitative agreement with the experiments 
on onset of superfluidity in 4 He in porous glasses, and on 
superconductor-insulator transition in homogeneous metallic 
films. 

The phenomenon of superconductor-insulator (SI) 
transition occurs in plethora of low-dimensional elec- 
tronic systems, examples ranging from Josephson junc- 
tion arrays M and homogeneous thin films §, D to high- 
temperature superconductors jl, and is believed to rep- 
resent a prototypical quantum (zero-temperature) phase 
transition. At low temperatures, as some controlling pa- 
rameter like thickness of a film is varied, the resistivity 
changes from a sharply decreasing to a continuously in- 
creasing function of temperature || . The good collapse 
of the resistivity data under scaling and near-universality 
of the critical value of the conductivity indicate a quan- 
tum (T = 0) critical point that separates two many-body 
ground states with different symmetries. A natural ques- 
tion arises: what is the mechanism of destruction of the 
superfluid ground state in a disordered system? One pos- 
sibly universal answer is provided by the so-called dirty- 
boson theory, which postulates that it is the loss of the 
superconducting phase coherence due to localization of 
Cooper pairs which is ultimately responsible for the SI 
transition [0], Q|, ||. Since on the scale of the diverg- 
ing phase-coherence correlation length Cooper pairs will 
appear as point particles, the SI transition would, under 
this hypothesis, in general fall into the same universal- 
ity class as the onset of superfluidity in A He in disor- 
dered media || , corrected for the long-rangeness of the 
Coulomb interaction. In principle, a way to assess the 
validity of this idea is to compare the measured criti- 
cal exponents with the calculations for the dirty-boson 
Hamiltonian. A strongly-coupled nature of the dirty- 
boson critical point, however, poses a fundamental obsta- 
cle to this procedure, and makes any but most qualita- 
tive understanding of the superfluid-Bose glass transition 
very difficult. The absence of a useful non- interacting 
starting point for disordered bosons forces one to rely on 
uncontrollable approximation schemes or turn to numeri- 
cal calculations [|l0| . This seems to be a common problem 
for the theories of interacting disordered low-dimensional 



quantum systems, apparent also in the fermionic systems 
of this type jy]. In fact, the dirty-boson Hamiltonian 
may be the simplest quantum problem that irreducibly 
contains the physics of interactions and disorder JTlfl , and 
as such has received a lot of attention through the years 

Recently, a new approach to the dirty-boson criticality 
has been suggested [|l3| , according to which the strongly 
coupled superfluid-Bose glass critical point in two dimen- 
sions (d = 2) could be understood as smoothly evolving 
from the zero-disorder critical point in d = 1 . The idea is 
to note that, by preventing the clean superfluid ground 
state in d = 1 to exhibit a true long-range order, the 
Mermin- Wagner theorem jl4j forces the SI fixed point in 
d = 1 to lie precisely at zero disorder pl| . In d = 1 + e a 
true long-range order becomes possible and the superfluid 
thus becomes more resilient to disorder, which causes the 
SI fixed point to shift to a finite, but small, value of 
disorder, controlled by the parameter e p3|. Although 
the dirty-boson transition probably lacks the upper crit- 
ical dimension [^6|, it has the lower critical dimension 
di = 1, and this in principle allows one to compute the 
universal quantities at the transition perturbatively in 
small parameter e = d — 1. Within this formalism, a 
particular symmetry of the low-energy action present in 
d = 1 guarantees that the dynamical critical exponent is 
z = d (z = 1 for Coulomb interaction) exactly ||13|| , in 
agreement with the expectation based on the compress- 
ible nature of the Bose-glass |l(| . The second, correlation 
length exponent v, however, turns out to be a perturba- 
tion series in yfk. On the experimental side, a directly 
measurable quantity is typically the product of the two 
exponents zv H , Q , and a meaningful comparison with 
experiment requires knowledge of the exponent v to some 
accuracy. In this Letter I present a field-theoretic method 
for higher-order calculations within the e-expansion for 
the superfluid-Bose glass transition, and use it to com- 
pute the correlation length exponent to two lowest orders 
in y/e. The result for both short-range and Coulomb in- 
teraction between bosons (see Table I) leads to values of 
v in d = 2 in a very good agreement with the experiments 
on the onset of superfluidity in 4 He in aerogel []9| and on 
the SI transition in thin metallic films Q , as well as with 
the Monte Carlo calculations on the dirty-boson Hamil- 
tonian [[l7|]. My calculation supports the idea that the 
SI transition in disordered electronic systems falls into 
the dirty-boson universality class, and establishes a way 
for a quantitative understanding of the SI criticality, as 
presently exists for the thermal critical phenomena p8| . 
The effort involved in higher-order calculations of the 
correlation length exponent and of the universal critical 
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Monte Carlo [|l7| 


e-expansion 




0.80 ±0.04 p| 


0.90 ±0.10 


0.81 


Zs 
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2.0 ±0.1 
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Vc 


1.2 ±0.2 24[ 


0.9 ±0.15 


1.03 


Z c 


1.0 ±0.1 § 


1.0 
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TABLE I. Comparison between experiment, Monte Carlo, 
and second-order e-expansion results for the critical ex- 
ponents, for short-range (s) and Coulomb (c) universality 
classes. 



conductivity is discussed. 

To be specific, consider the effective low- energy T = 
action for the disordered superfluid in d = 1 fll9f , | fL6f : 

K N f°° 

S=—J2 dxdT{c 2 [d x 6i(x,T)] 2 + [d T 9i{x,T)] 2 } (1) 
N ^,00 

-D^2 dxdTdT , cos2[6 i (x,T)-6 j (x,T')]. 
1,3=1 J -°° 

The coupling K is inversely proportional to the super- 
fluid density at some microscopic cutoff length A -1 , c 
is the velocity of low-energy phonons, and D is propor- 
tional to the width of the Gaussian random potential. 
The interaction between particles in (1) is taken to be 
short-ranged, and the standard limit on the number of 
replicas N — > is assumed. The field 8(x, r) is the dual 
phase Jll|, JhJ, and the above theory describes the de- 
struction of the second sound mode in the superfluid due 
to unbinding of topological defects (phase slips) at the 
point of transition in d = 1. Its role is similar to that of 
the sine-Gordon theory for the Kosterlitz-Thouless tran- 
sition in the 2D XY model |Q . To determine the macro- 
scopic state of the system, one is in principle interested 
in fate of the couplings K, c, and D as the cutoff in the 
theory is lowered. Under a change A — > A/s the com- 
binations of the coupling constants u = 3 — r/^ 1 , (where 
7] = Kc), k = 1/Kc 2 7 and W ~ D (to be precisely de- 
fined shortly), in d = 1 ± e dimensions are expected to 
renormalize according to the /3-functions pJ: 

u = e(u-3) + W + auW + 0(W 2 ), (2) 
W = uW + bW 2 + 0(uW 2 ), (3) 

k = (d — z)n, (4) 

where x = dx/dhi{s), and a and b are numerical coeffi- 
cients. The (i-dependent terms in the recursion relations 
(2)- (4) can be understood as deriving from the scaling 
of the superfluid density, p s f ~ K" 1 ~ £ 2_z_d , an j 
the compressibility, k ~ £ z ~ d |Q, pr| , where ^ is the 
diverging correlation length near the critical point, and 



z the dynamical exponent. Adopting the logic of the 
minimal subtraction scheme jl8| , the disorder-dependent 
terms in the recursion relations should be computed pre- 
cisely in d — 1, where one has the dual representation 
(1) of the low-energy theory on disposal. The symme- 
try of the interaction term in (1) under a transformation 
0i(x, t) — > 6i(x, t) ± f(x) for arbitrary f(x) implies then 
that there could be no disorder-dependent terms in Eq. 
(4) |l3|], so z — d at the fixed point. The correlation 
length exponent v follows from the linearization of the 
first two equations near the critical point at W* ~ u* ~ e. 
It is then straightforward to check that to the second or- 
der in e 1 / 2 

v, 1 = 3>e>+ K - '-e + 0{ei), (5) 

for short-range interactions. The 0(e 3 / 2 ) term follows 
from the higher-order terms in Eqs. (2)-(3). 

While the outlined procedure is conceptually simple, 
its implementation is made difficult by the fact that the 
action (1) has a compact form only in real space, and the 
interaction term contains all powers of the dual field. A 
similar obstacle appears in the calculation of the Koster- 
litz recursion relations beyond the lowest order in fu- 
gacity |^p[ . Here I will introduce a field-theoretic ap- 
proach to the problem, which also enables one to avoid 
the usual pitfalls of the momentum-shell renormalization 
group when applied beyond the lowest order. The gist of 
the method is to recognize that in d = 1, right at u = 
the coupling constant W becomes dimensionless, and the 
theory (1) appears to be just renormalizable. One then 
expects that the logarithmic divergences at the renormal- 
izable point d = 1 and u = will at small finite u show 
as poles when u — > 0. Since the coupling W acquires a 
finite scaling dimension for u ^ 0, the coefficients in the 
/3-functions are expected to stay finite as u — > 0. This 
is analogous to the standard procedure of dimensional 
rcgularization, commonly used to study thermal critical 
phenomena near the upper critical dimension Jis] , ex- 
cept that here the coupling constant u plays the role of 
dimension, while the real physical dimension is at first 
fixed at d — 1. When finally d —> 1 + e, the /3-functions 
are deformed into Eqs. (2)- (4). 

With this strategy in mind, consider the self-energy 
defined by the propagator of the dual phase in d = 1 
as G-\k,uj) = {2K/n)(u 2 + c 2 k 2 ) + £(w). It will prove 
useful to separate the first and the second order contribu- 
tions to S by writing it as S(w) = Si(cj)±E 2 (^)±0(-D 3 ), 
where £„ ~ D n . Simple calculation then gives 

/oo 
dr(l - e^Je-'M, (6) 
-00 

where the two-point correlation function /(t) = 
4(0(0, 0)(0(0, 0) - 0(0, t))) is given by the integral 

/• cA|r| dx 

f(r) = (3-u) -(1-e-*). (7) 

Jo x 



2 



When u — > 0, it readily follows that at small frequencies 
,2W ,1 



W' = W r + Su r W r + aW 2 



(15) 



E x (o;)= W 2 (- + 0(1)), 

7TC U 



(8) 



where I introduced the frequency-dependent, dimension- 
less coupling W = (4nD/(c 2 A 3 ))(cA/uj) u . After a te- 
dious but straightforward algebra one similarly finds 



71 ' El(w) ^ + /( W ), 



(9) 



where 



(10) 



I(oj) = 8D 2 / dydrdr' dv{\ - e luJT )e 

J OO 

[J?(y,t, J r,T / )(l + e- fc "'i(l-e- w )-l]. 
and f denotes a four-point correlation function: 

Wj T) t a = e -4((e(o,o)-fl(o ) T))(e( l /,«)-e( v ,«+T'))> _ ( n ) 

When u; — > 0, after rescaling the imaginary times and 
the length in the integral as lot —> r and uy/c — > y, 
the leading divergence in I (to) as u — + comes from the 
integration over small values of r and r'. To obtain the 
leading divergence in I(uj) it therefore suffices to expand 
F to the lowest order in r and r', to find that at small 
frequencies 



I{uj) = -uj 2 —W 2 {\ + 0(-)). 

ire u z u 



(12) 



The last equation is the central result of this work. Col- 
lecting all the terms back into the self-energy one recog- 
nizes the renormalized coupling rj r as the coefficient of 
<x> 2 -term in the propagator. In general, 



W W 2 ^,W 2 S 
Vr = V + — +x—+0( , 



(13) 



where the terms finite when u — > have been discarded, 
and x is a number determined from Eq. (12). After 
judiciously defining a renormalized disorder in (13) as 
W r = W + 2xW 2 /u, and rescaling it to bring the coeffi- 
cient of 0(W)-teim in the Eq. (2) to unity as 9W r — > W r , 
a differentiation with respect to ln(cA/w) leads to the 
Eqs. (2)-(3) for the renormalized couplings u r and W r , 
with the coefficients a = —2/3 and b = 2x/9, with 6 = 0. 
The subleading ~ W 2 /u term in the Eq. (13) determines 
the next, 0(W 2 ), term in (2), and the next-order correc- 
tion to v^ 1 . 

A remarkable feature of the perturbation series for v s 
is its independence of the renormalization procedure, i. e. 
of the non-universal finite parts of the self-energy which 
have been dropped in the last equation. To see this to 
the order of my calculation consider the most general 
redefinition of the coupling constants to the order W 2 



+ aW r + l3u r W r + ~fW 2 , 



(14) 



where the coefficients {a, ...ct} are finite, and dependable 
on the finite parts of the self-energy. It is easy to check 
that the recursion relations for the new couplings have 
the same form as the Eqs. (2)-(3), but with the coeffi- 
cients a' — a+a — S and b' = b — a+5. Interestingly, while 
the coefficients a and b by themselves are non-universal, 
the critical exponent depends only on their sum, which 
is perfectly universal, i. e. scheme independent. An in- 
variant similar to a + b appears also in the /3-functions 
for the sine-Gordon model p0| , where it determines the 
first correction to Kosterlitz-Thouless scaling. 

I expect the presented e-expansion to lead to divergent 
series; the point D = in the action (1) is non-analytic, 
since for D < the Gaussian probability distribution for 
the random potential becomes unbounded. Nevertheless, 
the hope is that the series in Eq. (5) for example, will be 
asymptotic, and that the few lowest terms may already 
lead to useful results. Indeed, estimating v s for d = 2 
from the simple sum of the first two terms gives v s = 0.81, 
within bounds found in the Monte Carlo calculations [0 . 
The experimental data of Crowell et al. on the onset of 
superfluidity and on the specific heat of He in aerogel at 
low temperatures H on its face value are consistent with 
the effective dimensionality of d = 2. Under this assump- 
tion the product of the two exponents in their experiment 
is zv — 1.60 ± 0.08. Assuming further that z s = 2 at the 
transition in d = 2 gives v — 0.80 ± 0.04. Although the 
uncertainty cited here should be taken with some reserva- 
tion, and the accuracy of the measurement falls short of 
the standards in thermal critical phenomena, the result 
appears to be in excellent agreement with my calculation 
(see Table I). It is worth noting that the inequality pi] 
v>2/d seems to be violated both by the experiment and 
by my estimate. It has been argued recently |^2j that the 
above inequality is an artifact of the particular averag- 
ing procedure, and that the true exponent is in fact not 
bound from below. It would be interesting to see if the 
higher-order corrections eventually push the value of v s 
above unity in d — 2, or indeed v s < 1 as the experiment 
and the present calculation suggest. 

To make a comparison with the experiments on SI 
transition in homogeneous thin films with thickness as 
the tuning parameter 0] it is necessary to take into ac- 
count the long-range Coulomb interaction between the 
Cooper pairs. As explained in detail elsewhere [^3), |23] ] 
within the present scheme this may be simply accom- 
plished by defining the Coulomb interaction as V c (r) = 
e 2 J d d q exp(iq-f)/ g rf_1 , so to coincide with the (5-function 
repulsion in d = 1. The only change in the calculation 
then is that z c — 1, and that e — > e/2 in the Eq. (2), while 
the disorder-dependent terms in the recursion relations, 
which follow from d = 1, remain the same. The first two 
terms in the series then give v c — 1.03, in accord with the 
Monte Carlo results |l7j (see Table I). Experiment finds 
z c = 1.0±0.1, by suppressing the transition temperature 
with the magnetic field or by scaling of resistance with 
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the electric field ||. Collapsing the resistivity data ]24| ] 
then gives the experimental value of v c = 1.2 ±0.2, again 
in a very good agreement with my result. 

As mentioned earlier, the next term in the series (5) re- 
quires only the computation of the subleading, 0(W 2 /u) 
term in Eq. (13). Here, however, it appears that it is no 
more enough to know the correlator F (after rescaling 
the lengths with u>) only at small rescaled r and r', as 
it was for the leading divergence in Eq. (12). In light of 
the likely asymptotic nature of the expansion, this is left 
for future work. 

Another universal quantity of interest is the critical 
conductivity in d — 2 || , which, aside from the universal 
unit el/h, for bosons of charge e* can be obtained as a 
Laurent series in e The lowest order term was ob- 

tained in ^3|, and for d — 2 the result <j c = 0.69e^/h 
agrees with the low-temperature experiment on Bismuth 
films |24|] quite well. It would nevertheless be useful to 
compute the next-order correction, and see if it pushes 
the result towards the self-dual value of el/h, to which a 
large number of experimental results seems to converge. 
Calculating the next-order term in the critical conduc- 
tivity would in principle proceed along the same lines as 
here, except that one needs to perform the analytic con- 
tinuation to real frequencies to obtain the real-time dc 
res pon se at low temperature. This problem was solved 
in p3j to the lowest order, but applying the same trick 
in the next order term seems not straightforward. A pre- 
liminary analysis also suggests that the second-order con- 
tribution to the universal conductivity requires a compu- 
tation of both leading and subleading divergences in the 
second-order contribution to the self-energy. 

In conclusion, it is shown that the expansion around 
the lower critical dimension di — 1 for the superfluid- 
Bose glass critical point allows a field-theoretic formula- 
tion that facilitates a systematic higher-order calculation 
of the critical exponents. The computation to two low- 
est orders yields results for the correlation length crit- 
ical exponents in d = 2 in respectable agreement with 
the experiment and Monte Carlo calculations, both for 
the short-range and the Coulomb interaction between 
particles. The results suggest that the superconductor- 
insulator transition in homogeneous thin films is in the 
universality class of disordered bosons. 
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